Role of a parallel magnetic field in two dimensional disordered clusters 
containing a few correlated electrons 
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An ensemble of 2d disordered clusters with a few electrons is studied as a function of the Coulomb 
energy to kinetic energy ratio r s . Between the Fermi system (small r s ) and the Wigner molecule 
(large r s ), an interaction induced derealization of the ground state takes place which is suppressed 
when the spins are aligned by a parallel magnetic field. Our results confirm the existence of an 
intermediate regime where the Wigner antiferromagnetism defavors the Stoner ferromagnetism and 
where the enhancement of the Lande g factor observed in dilute electron systems is reproduced. 

PACS: 71.10.Hf, 71.10.-w, 71.30,+h, 73.23.-b 



The magnetization of quantum dots is the subject of 
recent studies since their Ohmic resistances are mea- 
sured as a function of a gate voltage in the Coulomb 
blockade regime. The Pauli rule (electrons populating 
the orbital states of a dot in a sequence of spin up - spin 
down electrons) leads to a bimodal distribution of the 
conductance peak spacings. The possibility of a spon- 
taneous magnetization S of their ground state due to 
electron-electron interactions has been proposed to ex- 
plain the absence of such a distribution in the experi- 
ments. Since the Coulomb energy to Fermi energy ratios 
r s are not too large in those dots, the interactions are 
mainly described assuming the Hartree-Fock (HF) ap- 
proximation. A break-down of the Pauli rule is obtained 
if, by placing electrons at higher orbitals, the gain in ex- 
change energy exceeds the loss in kinetic energy. This is 
the Stoner mechanism || which eventually gives a ferro- 
magnetic ground state. The mesoscopic fluctuations of 
S are usually studied PJ assuming some random matrix 
toy model valid in the zero dimensional limit. We nu- 
merically study in this letter another limit where the 2d 
character of the dot is kept and where the interactions 
are exactly taken into account. Therefore, our study is 
directly relevant for quantum dots created from dilute 2d 
electron gases and may give insights into the physics of 
the dilute 2d electron gas (strongly correlated limit). 

The dilute limit is opening a new frontier for experi- 
ment and theory. 2d electron gases can be created in 
field effect devices at low carrier densities n s and large ra- 
tios r s oc ris^ 1 / 2 . The dependence of the 2d resistance R 
as a function of the temperature T and the applied mag- 
netic field B has been studied up to large values of r s . 
Wigner crystallization, which remains to be confirmed 
in the experiments, is seen in computer simulations for 
values of r s where a metal-insulator transition (MIT) is 
observed [^|,^| . One has an insulator when r s is large 
and a metal when r s is smaller. The MIT occurs |?]] 
at r s w 10 in disordered devices, at larger r s in cleaner 
devices. Since Anderson localization takes place when 



r s — > 0, the metallic behavior of R(T) observed for in- 
termediate values of r s suggests the possible existence of 
an intermediate metallic phase between two insulating 
limits of different nature: a Fermi system of Anderson 
localized states on one side and a pinned Wigner crys- 
tal on the other side. The existence of such a phase is 
supported by exact studies ||[| of small disordered clus- 
ters where three different regimes are observed for the 
persistent currents. Nevertheless it remains unclear that 
the observed intermediate regime is not a simple finite 
size effect, and that it does correspond to a new kind 
of metal. An analogy with the physics of 3 He has been 
recently discussed [jl(| . 

Renormalization group equations |p"lf obtained for 
weak disorder and weak interaction underline the role 
of the spins and suggest that a 2d metal might exist. 
The behavior of R when a parallel magnetic field B is 
applied (no orbital effects) confirms the role of the elec- 
tronic spins. For the intermediate r s where R(T) be- 
haves as in a metal, a positive magnetoresistance Q 
saturating at a large value for B > B sat is observed. 
Small-angle Shubnikhov-de-Haas measurements |T^ | con- 
firm that B sat is the necessary field to polarize the elec- 
trons completely. Following ref. [fl3||, B sat oc (n s — n c ) 
where n c is the density where the MIT takes place at 
B = 0. This behavior is observed for densities n s corre- 
sponding to 3 < r s < 10 and gives g(n s ) oc n s (n s — n^ 1 
for the Lande g factor. A shift of n c (defined from the 
I —V characteristics) is also observed Q : n c increases 
as a function of B before saturating when B > B sat . 
This typically gives for the critical ratios r c s : r c s {B = 
Q)/r c s {B > B sat ) » 1.2. 

Our study of an ensemble of small disordered clusters 
confirms that the intermediate regime observed j§) using 
spinless fermions remains when the spin degrees of free- 
dom are included. When the electrons are not polarized 
by a parallel field B, the ground state begins to be de- 
localized by the interaction before forming the localized 
Wigner molecule. When the spins are aligned by B, this 
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derealization effect is suppressed and the crystallization 
threshold is shifted by a small amount compatible with 



the experimental observations 14 1. The polarization field 



B sat cx r~ 2 found in our simulations is consistent with 
the law B sat cx n s found in Ref. |l3| in a similar interval 
of intermediate values for r s . Magnetic signatures of the 
intermediate regime are given. The local magnetic mo- 
ments (S — 1) appearing above the mesoscopic Stoner 
threshold (r s ss 0.35 for the studied clusters) begin to 



disappear above a first threshold 
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2.2 where they 



have a maximum probability. Above rf s , the Stoner 
mechanism and hence the HF approximation break down, 
B sat tx rj 2 and the field B necessary to create a local mo- 
ment {S — 1) in a S = cluster becomes independent of 
r s . Those behaviors disappear around rY s ~ 10 where 
the electrons crystallize and form an antiferromagnetic 
(S = 0) Wigner molecule. 

The clusters are described by an Hamiltonian H of N 
electrons free to move on a L x L square lattice with on- 
site disorder and periodic boundary conditions (BCs). 
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The operators Ci t „ (cj a ) destroy (create) an electron of 
spin g at the site i and rii^ — c\ a Ci. a . The hopping term 
—t couples nearest-neighbor sites and the random po- 
tentials Vi are uniformly distributed inside [— W/2, W/2] 
with W — 5 (diffusive regime for non interacting parti- 
cles). Two electrons on the same site costs a 2U Hubbard 
energy and a U/\i — i'\ Coulomb energy if they are on two 
sites i and i' separated by \i — i'\ (shortest length separat- 
ing i and i' on a square lattice with periodic BCs). The 
use of the dimcnsionless ratio r s = U/(2ty/nn s ) where 
n s = N/L 2 allows us to compare our results obtained as 
a function of U and the experimental results obtained as 
a function of n s . 

We denote 5 and iS 2 the total spin and its compo- 
nent along an arbitrary direction z. Since [<S 2 ,7Y] = 
[iS z ,7i] = 0, Ti can be written in a block-diagonal form, 
with N + 1 blocks where S z = -N/2, N/2 respec- 
tively. The blocks with S z and — S z are identical in the 
absence of a magnetic field B. When 5 = 0, there is 
no preferential direction and the groundstate energy Eq 
does not depend on S z . For a groundstate of total spin 
S, H has 2S+ 1 blocks with the same lowest eigenenergy 
E a (S 2 ) since E a (S 2 ) = E {S 2 ,S Z ); s,=-S,-S+i,...,S-i,s. 
Therefore, the number Nb of blocks of different S z and of 
same lowest energy gives the total spin S — (Nb — l)/2 
of the groundstate. We consider N — 4 and L = 6, as 
in ref. [11. S z — 2,1,0,-1,-2 give two characteristic 
energy differences Ai = E (S Z = 1) — E (S Z = 0) and 
A 2 = E (S Z = 2) - E (S Z = 0). The E (S Z ) energies are 



exactly obtained using Lanczos algorithm. The size of 
the matrices to diagonalize are N(S Z ) = 396900,257040 
and 58905 for S z = 0, 1, 2 respectively. 

A parallel magnetic field B defines the z-direction 
and removes the S z degeneracy by the Zeeman energy 
—gnBS z . The ground state energy and its magnetiza- 
tion are then given by the minimum of Eq(S 2 , S z , B = 
0) — g/iBS z . For a S = groundstate without field, the 
value B* for which E Q (S Z ) - g^B*S z = E (S Z = 0) de- 
fines the field necessary to polarize the system to S > S z . 
The total A 2 and partial Ai polarization energies are the 
Zeeman energies necessary to yield 5 = 2 and S = 1 re- 
spectively for a cluster with S = 0. 

When r s = 0, the two one body states of lowest en- 
ergy are doubly occupied and S = (S — 1/2 if N is 
odd). To polarize the S = ground state to S = 1 cor- 
responds to the transition of one electron at the Fermi 
energy and costs an energy equal to the one body level 
spacing. p(Ai) is then given by the spacing distribution 
p(S) between consecutive one body levels, the Wigner 
surmise (nS/2) cxp — (ttS 2 /A) in the diffusive regime. A 2 
corresponds to the sum of a few one body excitations. 

When r s is large, the 4 electrons occupy the four sites 
c(j) j = 1, . . . , 4 of the square configuration of side a = 3 
with the lowest substrate energy J2^=i v c(j)- 11 |0 > de- 
notes the vacuum state, the ground state in this limit 
becomes \^! c >= Ylj=i c I(j) a \® > with a spin indepen- 
dent energy E c . This square can support 2 N — 16 spin 
configurations. We summarize the main results of a per- 
turbative expansion of E c in powers oft/U. The spin de- 
generacy of E c is removed by terms of order t(t/[/) 2a_1 , 
which is the smallest order where the 16 spin configu- 
rations can be coupled via intermediate configurations 
allowing a double occupancy of the same site. There- 
fore, 2a — 1 is the order where the perturbation begins 
to depend on S z and Ai as A 2 cx t{t/\J) 2a ^ 1 — > when 
t/U — > (we have numerically checked this decay when 
r s > 100). Moreover, the correction to E c depending 
on S z and oc t(t/U) 2a ~ 1 is given by an effective antifer- 
romagnetic Heisenberg Hamiltonian. The S = ground 
state for large r s correspond to 4 electrons forming an an- 
tiferromagnetic square Wigner molecule. This does not 
reproduce the ferromagnetic ground state suggested in 
Ref. jn|. However A 2 is very small when r s is large, 
and the antiferromagnetic behavior can be an artefact 
due to the square lattice. Without impurities, a quasi- 
classical WKB expansion []l6| shows that 3 particle ex- 
changes dominate in the continuous limit, leading to fer- 
romagnetism. Recent Monte-Carlo calculations |l7]] sug- 
gest that the crystal becomes a frustrated antifcrromag- 
net closer to the melting point. 

The perturbative corrections oc t(t/U) 2a ~ 1 depend on 
the random variables via 

11 j% (Ec-Ej)- 1 where the 
Ej are the classical energies of the intermediate configu- 
rations. Ej is the sum of an electrostatic energy and of 
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a substrate energy E S (J) = J2k=i v J(k)- Due to the high 
order 2a— 1 of the correction, a normal distribution for Ej 
leads to a log- normal distribution for YlJ=i {-® c — Ej)~ 1 . 
Therefore p(Ai) and p(A2) should be log-normal when 
r s is large. 



1.0 - 



0.5 - 



0.0 



0.3 - 



0.0 



p 1 I 1 I 1 

i P(5,) 




*, -log(A) 
















i 


K i 

VSe i i • • • • •! 


,\g — V« t . i 




JJW. 



- 1.5 

1.0 
P 

0.5 



- 0.6 



0.0 



2 3 4 5 

-log(A) 



3 4 

-log(A) 



o.o 



FIG. 1. Distributions of Ai and A2 at different values 
of r s . Upper left: p(Si) at r s = (circle) and 2.5 (square) 
where Si = Ai/ <Ai>. The continuous line is the Wigner 
surmise. p(— logAi) (circle) and p(— logA2) (diamond) at 
r s — 2.5 (upper right) 5.8 (lower left) and 16.8 (lower right) 
respectively. The thick bars (put at right edge of the figures) 
give the peaks <S(Ai) of the bimodal p(Ai). The continuous 
lines are normal fits. 

p(Ai) is given in Fig. |l] for different r s . The expected 
Wigner surmise takes place for r s = 0. A small interac- 
tion quickly drives p(Ai) towards a bimodal distribution, 
with a delta peak at Ai = and a main peak centered 
around a non zero value of Ai. The delta peak gives 
the probability to have spontaneously magnetized clus- 
ters with 5 = 1. The main peak gives the field B nec- 
essary to create 5 = 1 in a cluster with 5 = 0. The 
logarithmic scale used in Fig. [l] underlines the bimodal 
character of the distribution and confirms that the main 
peak becomes log-normal when r s is large. The distri- 
bution of A2 is not bimodal: a fully polarized cluster 
has never been seen when B = 0. A2 becomes also log- 
normally distributed when r s is large. 

In Fig. ||, the fraction M of clusters with 5 = 1 at 
B = is given as a function of r s . One can see the meso- 
scopic Stoner instability || taking place at r s w 0.35. 
The Stoner mechanism should eventually give fully po- 
larized electrons. This is not the case, the increase of 
M breaks down when r s = rf s « 2.2, a value where 
the Stoner mechanism and hence the HF approximation 
break down. In the same clusters, the HF approximation 
fails p8|] to describe the persistent currents of 4 spin- 



less fermions when r s > rf p w 5. rf p takes a smaller 
value rf s when the spin degrees of freedom are included. 
Above rf s , M regularly decreases to reach a zero value 
for r^ s w 9 where an antiferromagnetic square molecule 
is formed. In the intermediate regime, there is a competi- 
tion between the Stoner ferromagnetism and the Wigner 
antifcrromagnetism. Since the 5 = clusters are char- 
acterized by log-normal distributions, the ensemble av- 
erages < logAi > and < logA 2 > (without taking into 
account the 5=1 spontaneously magnetized clusters) 
define the typical fields B necessary to yield 5 = 1 or 
5 = 2ina5 = cluster. Fig. ^| provides two magnetic 
signatures confirming the existence of a novel intermedi- 
ate regime between the Fermi glass (r s < rf ) and the 
Wigner glass (r s > rf s ): < logAi > becomes roughly 
independent of r s , while A2 cx r~ 2 . Very remarkably, this 
rj 2 dependence is consistent with the n s dependence of 
B sat seen in the experiments [O . 




FIG. 2. As a function of r 3 , fraction M (diamond) of 
clusters with S = 1 at B = 0, partial < log Ai > (circle) and 
total < log A2 > (triangle) energies required to polarize S — 
clusters to S = 1 and 5 = 2. The straight line corresponds to 
0.25 - 2 log r s . 

We now study three quantities calculated from the 
three eigenstates l^ol^z) > of lowest energy with 
5z = 2,1,0 respectively: (i) the number £(5 Z ) = 
N 2 {J2i Pi (Sz))^ 1 of occupied lattice sites (p l {S z ) =< 
*o(5^)| t2 a n i,cr\^o(Sz) >); (h) the crystallization pa- 
rameter 7(5 Z ) = max r C(r, S z ) — min r C(r, S z )) where 
C(r,S z ) = N- 1 Y Jl p l {S z ) Pl+r {S z ). i{S z ) = (1) if the 
state is a liquid (a crystal); (iii) the longitudinal Ii(S z ) 
and transverse It{S z ) components of the total persis- 
tent current driven by an Aharonov-Bohm flux <f> = n/2 
enclosing the 2d torus along the longitudinal direction 
(cf) = it corresponds to anti-periodic longitudinal BCs 
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and periodic transverse BCs, see Ref. ||). 

£(S Z ) (shown for a single cluster in Fig. ||) depends 
on S z for small r s and becomes independent of S z for 
large r s . At r s = 0, $o(S z = 2,1,0) occupy respec- 
tively 4, 3, 2 one body states while the Wigner molecule 
occupies 4 sites only at large r s . The ensemble average 
< £,(S Z — 2) > is maximum when r s — and decays as 
r 8 increases, suggesting the absence of derealization for 
the polarized system. The non polarized system behaves 
differently, since < £(S Z — 1,0) > first increase to reach 
a maximum «< £(S Z = 2, r s = 0) > before decreasing. 

10 20 30 




FIG. 3. As a function of r 3 , numbers £(S Z ) of occu- 
pied sites for a typical sample (upper left); ensemble aver- 
ages < £(S Z ) > (upper right) and crystallization parameter 

< y(S z ) > (lower left) for the lower energy states with S z — 
(circle), S z — 1 (square) and S z = 2 (triangle); Average lon- 
gitudinal < Ii(S z = 0) > (circle) and transverse (diamond) 

< h{S z = 0) > components of the total persistent current 
(lower right). 

In the lower left part of Fig.|| one can see that charge 
crystallization is easier when the clusters are polarized 
than otherwise. The shift of the critical threshold rY is 
consistent with the shift of the critical density reported 
in Ref. [Q. The arrows indicated in Fig. || underline two 
consequences of a parallel field B: smaller number of oc- 
cupied sites and smaller crystallization threshold. This is 
qualitatively consistent with the large magnetoresistance 
observed in the metallic phase for intermediate r s . Fig. 
H (lower right) gives the total longitudinal < h(r s ) > 
and transverse < It{r s ) > currents for S z = 0. When the 
spin degrees of freedom are included, the conclusions pre- 
viously obtained for spinless fermions j^^] remain valid: 
there are a Fermi regime where longitudinal and trans- 
verse currents coexist, an intermediate regime where the 
transverse current is suppressed while the longitudinal 



current persists, and a Wigner regime with vanishing per- 
sistent currents p9| . 

In summary, our study suggests that Coulomb repul- 
sion by itself provides a qualitative mechanism able to 
give together a derealization of the electrons ( "metallic 
behavior") and a large positive magneto-resistance for 
intermediate value of r s . The shift of the critical thresh- 
old as a function of B and the dependence of B sat as 
a function of r s are reproduced. Two important ques- 
tions remain opened. Do the finite size behaviors which 
we have observed persist in the thermodynamic limit? 
What are the physical explanations of the interplay be- 
tween the Stoner ferromagnetism and the Wigner antifer- 
romagnetism seen for intermediate values of r s . An inter- 
mediate regime where a floppy antiferromagnetic Wigner 
molecule with vacancies coexists |2(| with a ferromag- 
netic liquid of intersticial particles is suggested as a pos- 
sible explanation. 
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